Orbital Angular Momentum Origin of Rashba-type Surface Band Splitting 
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We propose that existence of local orbital angular momentum (OAM) on the surfaces of high- 
Z materials plays a crucial role in the formation of Rashba-type surface band splitting. Local 
OAM state in a Bloch wave function produces asymmetric charge distribution (electric dipole). 
Surface-normal electric field then aligns the electric dipole and results in chiral OAM states and 
the relevant Rashba-type splitting. Therefore, the band splitting originates from electric dipole 
interaction, not from the relativistic Zeeman splitting as proposed in the original Rashba picture. 
The characteristic spin chiral structure of Rashba states is formed through the spin-orbit coupling 
and thus is a secondary effect to the chiral OAM. Results from first principles calculations on a 
single Bi layer under an external electric field verify the key predictions of the new model, including 
the direction of the spin texture which we predict to be in opposite direction to the conventional 
Rashba picture. 
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In the periodic band of a solid with inversion and time- 
reversal symmetries, Kramer's theorem mandates that 
each momentum state at k be spin degenerate [1]. The 
spin degeneracy may be lifted, however, on surfaces of 
solids or interfaces of hetero-structures where the inver- 
sion symmetry is broken [2]. Lifting of the spin degener- 
acy due to inversion symmetry breaking (ISB) is usually 
called the Rashba effect [3], which also entails a chiral 
spin structure along a constant energy contour as a con- 
sequence of the spin-momentum locking. Surface energy 
splitting and concomitant chiral spin structure have been 
experimentally observed on Au(lll)[4-7], Bi[8], Sb[9], 
and on some alloys as well[10]. The effect has drawn 
additional attention recently due to its relevance for the 
surface states of topological insulators [11]. 

In spite of its phenomenological success, the Rashba 
Hamiltonian Hr = Ar<t • (kx z), Ar= Rashba parameter, 
<x=Pauli matrix, k=electron momentum and i=surface 
normal, has several shortcomings when applied to solid 
surface phenomena. Estimates of Rashba energy Er us- 
ing the realistic work function at the surface gives Er ~ 
10~ 6 eV[12], a value far too small to account for the ob- 
served energy splitting as large as a few hundred meV[4- 
9] . In addition, the Rashba picture has difficulty explain- 
ing the observation that the splitting increases with the 
atomic spin-orbit coupling (SOC) strength a[5, 6] be- 
cause the strength of the surface electric field should not 
vary too much as a function of the atomic number [13]. 

Several recent studies have presented improvements 
and/or alternatives to the original Rashba model. Within 
the tight-binding model, Petersen et al. showed that the 
split energy indeed should be proportional, simultane- 
ously, to a and surface potential gradient 7 [12]. Even 
though not explicitly discussed in their work, the energy 



splitting also comes out to be proportional to the electron 
momentum k, which is essential in forming a Dirac cone- 
like dispersions. Some recent first-principles calculations 
show that the energy splitting is closely related to the 
electric potential energy of the asymmetric charge distri- 
bution at the surface[14, 15]. Orbital-mixing character 
of surface bands was emphasized by Bihlmayer et al[16\. 
Surface band splitting can be very large purely by band 
structure effect at specific wave vectors [17]. However, ef- 
forts to find the microscopic origin of the surface energy 
splitting and its consequences so far have been largely 
numerical, and do not seem to confront the underlying 
physical mechanism. The tight-binding model analysis in 
Ref. 12 points out the simultaneous presence of SOC and 
surface electric field as key factors in the Rashba split- 
ting, but does not incorporate the wave function asym- 
metry found in later band calculations. The physical ori- 
gin of the chiral spin structure was not addressed in band 
calculation approaches[18], and one needs to resort to the 
original Rashba picture when the chiral spin structure is 
discussed. In addition, as we will discuss later, the direc- 
tion of the spin texture in strong SOC cases predicted by 
the Rashba picture is opposite to what has been experi- 
mentally measured, which has been somehow unnoticed. 

In this article, we argue that the most important aspect 
of the energy splitting and resulting chiral spin structure 
is the existence of orbital angular momentum (OAM) in 
high-Z materials due to the strong atomic SOC. We fur- 
ther show that such an OAM state in combination with 
the electron momentum causes asymmetric charge dis- 
tribution, which in turn determines not only the energy 
level but also the OAM direction through electric dipole 
interaction. The spin angular momentum (SAM) direc- 
tion is determined by the OAM direction because of the 
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FIG. 1: Energy levels of p states (a) in an atom, (b) with a 
weak crystal field, (c) with a strong SOC and (d) in a tight 
bing state. Hcf and E s are the crystal field Hamiltonian and 
surface electric field, respectively. 



strong SOC, forming a chiral SAM structure. Therefore, 
SAM texture is a secondary effect due to the OAM chi- 
ral structure, not the primary effect, contrary to what 
has been been viewed so far. In order to verify the new 
model, we also present first-principles calculation results 
on a single-layer Bi sheet under an external electric field. 
The results are consistent with the new model. 

We consider a system consisting of p-orbitals in this pa- 
per, which is of practical importance as the relevant sur- 
face bands of elements such as Bi, Sb, and Pb exhibiting 
Rashba-type splitting consist mainly of p-orbitals. Even 
the valence states of Au surface show p-orbital character 
due to strong mixing between s- and p-orbitals by surface 
electric field[14]. 

For an intuitive understanding, we consider a case in 
which the SOC is very large. Energy levels for such case 
are schematically shown in Fig. 1. Without SOC, OAM 
is quenched in the presence of crystal-field and, p-orbitals 
are the energy eigen-states as illustrated in the figure. 
When SOC is turned on, the total angular momentum J 
eigen-states become the energy eigen-states instead. The 
most important aspect of this is that OAM comes back 
to the picture as a relevant low-energy degree of free- 
dom as will be explained later. At the last stage of Fig. 
1, the OAM in combination with the electron momen- 
tum k causes asymmetric charge distribution or electric 
dipole moment whose direction is opposite for the states 
with opposite OAM directions. In the ISB circumstances 
such as on surfaces, the electric field splits the states 
through electric dipole interaction. Subsequently, up- 
per and lower Dirac cone (UDC and LDC, respectively) 
states depicted in the figure are formed. 

The first step to understanding the process is to see 
how the asymmetric charge distribution is formed. Let's 
focus on the J = 1/2 doublet 
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FIG. 2: (a) A tight binding state built from J1/2 states. Each 
J1/2 state is the local state at each atom. Phase of each state 
is represented by the sinusoidal grey line. Electron density 
of a tight binding state integrated along x-direction for (b) 
k = O.Itt & n = x, (c) k = 0.05tt &h = x, and (d) k = O.Itt 
& n = — x. The ratio between Bohr radius for the p-orbital 
and the lattice constant is taken to be 7/20. Black dots in 
(b)-(d) indicate positions of atoms, (e) Integrated (along both 
the x- and y-directions) electron density as a function of z for 
the three cases. 
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which can be used to form a coherent state \n) satisfying 
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Here n = (sin 9 cos 0, sin 9 sin </>, cos 9) is the unit vec- 
tor indicating the direction of the total angular momen- 
tum average (J), J = L + (l/2)er. A Bloch eigenstate, 
characterized by both the wave vector k and the to- 
tal angular momentum orientation n, can be formed, 
|n, k) = N~ x / 2 Y^,j e lk ' ri |n,ri), as schematically shown 
in Fig. 2(a) for J (thus L too) parallel to the x-direction. 



Here, N is the number of sites and \n, Tj) refers to the 
Wannier state at the atomic site r^. One can show, as in 
Fig. 2, that the density distribution of the Bloch state 
depends on both k and h now, and in particular results 
in nonzero dipole moment that varies with them both. 

To demonstrate that OAM can induce nonzero dipole 
moments, we plot in Fig. 2(b)-(d) the density of the 
Bloch state for different k and fi values. For simplicity, 
we use hydrogenic 2p states as the orbital part of the 
Wannier state. The plotted densities are integrated elec- 
tron densities of the Bloch state along the x-direction. It 
is clear from Fig. 2(b) that electron density is higher 
in the — z region. This asymmetric electron distribu- 
tion can be understood in the following way. One can 
view the tight-binding state as a superposition of a free 
electron and local OAM states. When the phase veloci- 
ties of the electron and OAM states are pointing in the 
same direction (— z region), there is a constructive inter- 
ference which results in a higher electron density while 
opposite motions cause destructive interference, yielding 
lower electron densities (+z region). 

Electron densities for other combinations of k and n 
values are plotted in Figs. 2(c) and 2(d). The charge 
density becomes less asymmetric as the momentum k de- 
creases from 0.17T to 0.057T. It should eventually become 
symmetric at k=0 for an obvious reason. On the other 
hand, the density is higher in the +z region when angular 
momentum is flipped to n = -1 state (Fig. 2(d)). The 
trend can be seen more clearly by comparing integrated 
electron densities for the three cases plotted in Fig. 2(e). 
One can intuitively understand that the resulting electric 
dipole moment d should be proportional to L x k. 

The significance of the OAM induced electric dipole 
moment is that the surface electric field E s , directed nor- 
mal to the xy plane, couples to the dipole moment and re- 
sult in the splitting of energy among different L (thus n in 
equation (2)). To see this, we can compute dipole inter- 
action energy -Ed(k) = — eE s J dr z|V>7i,k(i")| 2 , e=electric 
charge, isoelectric field strength, and Vri,k( r )=wave 
function of the Bloch state, to find 



Ed(n, k) oc — Wz ■ (h x k) = Wk ■ (n x z) 



(3) 



in the small-k limit, and W — eEa denotes the work 
function of the surface. Our non-relativistic, electrostatic 
consideration leads to the correct form and the energy 
scale of the Rashba-type splitting. 

An important consequence of the dipole interaction 
is that L now has a preferred direction. The maxi- 
mum/minimum energy is obtained for L x z directed 
either parallel or anti-parallel to the k-vector (see Fig. 
3(a)). Therefore, the direction of the L is determined by 
the relationship E s || ±L x k, and it is the OAM that 
determines the locking of the angular momentum to the 
electron momentum. Due to the strong SOC, SAM S 
is anti-aligned to L, which is a secondary effect. The 
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FIG. 3: (a) Interaction of the electric dipole moment and 
surface electric field E s , and resulting alignment of the OAM. 
The lobes schematically show the electron distribution of the 
tight binding state. This particular configuration represent 
the lower energy state of Ji/2- (b) Resulting SAM and OAM 
textures as well as the (Dirac) band dispersion, (c) Angular 
momentum structures predicted by the Rashba model (left) 
and new model (right). 



resulting SAM and OAM structures for J = 1/2 are il- 
lustrated in Fig. 3(b) which are compatible with earlier 
reports[19]. Here, we point out that E s originates from 
the ISB, and therefore wave functions are determined by 
the ISB. 

There are two crucial differences between the SAM 
structures predicted by the two models. First of all, 
the SAM structure in the Rashba picture for strong 
SOC cases is found to be opposite to that in the new 
model shown in Fig. 3(b). For an electron moving in 
the jc-direction in Fig. 3(a), the effective magnetic field 
'B e ff = —{fi/mc 2 )k x E s is along the +y-direction. For 
such magnetic field, spin along the — y-direction gives the 
lower energy. This results in clock-wise chiral spins for 
the lower Dirac cone, opposite to that of the new model 
as shown in Fig. 3(c). Experimental data from, for exam- 
ple, Au(lll) surface states which correspond to J = 1/2 
states are consistent with the result of the new model. 
This problem in the Rashba picture has not been noticed 
so far. The other difference is that the SAM structures 
for the J = 1/2 and 3/2 cases are opposite in the new 
model. This is because OAM and SAM are directed par- 
allel to each other for the J = 3/2 states. On the other 
hand, the OAM structure is expected to be the same for 
all the Dirac cones (split bands near the V point) as the 
Rashba Hamiltonian does not consider the OAM. The an- 
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FIG. 4: (a) Calculated band struture along the high sym- 
metry directions for a single Bi layer in a triangular lattice 
under an external electric field of 3 V/A along the z direction 
as shown in the inset, (b) - (g) SAM and OAM expectation 
values. Bands are numbered in panel (a). Red and blue ar- 
rows represent SAM and OAM. The k-space region for (b)-(g) 
is shown as the dotted square in the first Brillouin zone on 
the right-hand side of panel (a). 



gular momentum structures predicted by the two models 
for an electric field along the z are schematically shown 
in Fig. 3(c). Confirmation of the SAM structures would 
be a definitive evidence for the new model. 

One may think about checking out the SAM struc- 
tures on the known surface states such as Au(lll) sur- 
face states. However, J = 3/2 states are normally in the 
bulk states and cannot be separated from the bulk states. 
In addition, there could be a question on the direction of 
the surface electric field even though it is more likely out- 
of-surface direction. In order to circumvent these prob- 
lems, we performed the first principles calculation on a 
single-layer of Bi in a triangular lattice under an external 
electric field. In this way, we simulate the surface with 
a definitive field direction and without bulk states (inset 
of Fig. 4(a)). The applied field (surface field in the z 
direction) is 3 V/A which is a reasonable number con- 
sidering the energy (work function) and length (^atomic 
size) scales. 

For the density-functional theory (DFT) calculations 
within the local-density approximation, we used the DFT 
code, OpenMX,[20] based on the linear-combination-of- 



pseudo-atomic-orbitals (LCPAO) method[21] and spin- 
orbit couplings were included via the norm-conserving, 
fully relativistic j-dependent pseudopotential scheme in 
the non-collinear DFT formalism[22-24]. LCPAO coeffi- 
cients at specific /c-points were used to calculate the SAM 
and OAM. 

The resulting band structure is plotted in Fig. 4(a). 
The initially degenerate bands are split upon application 
of the field as expected. It should be noted that the top 
bands (5 and 6) primarily consist of p x and p y orbitals 
and are less susceptible to formation of OAM. This re- 
sults in a much smaller splitting and thus supports our 
view that formation of OAM is essential in the energy 
splitting. On the other hand, split energies of other bands 
are similar to that of the Bi surface states even though a 
detailed level comparison is not meaningful. 

Calculated expectation values of SAM and OAM near 
the r point are plotted in Figs. 4(b) - 4(g). We first 
note that OAM in bands 5 and 6 (where the splitting 
is very small) is much smaller than that in other bands, 
which again supports the new model. In comparing the 
results with the structures in Fig. 3(c), one finds that 
1) the SAM direction for J — 1/2 in the Rashba picture 
is indeed wrong, 2) SAM directions for J = 1/2 and 
3/2 cases are opposite, suggesting that chiral structures 
are determined by the OAM, and 3) SAM and OAM 
are anti-parallel and parallel for J = 1/2 and 3/2 cases, 
respectively. These results are exactly what are expected 
in the new model, and thus conclusively prove our new 
understanding. 

What has not been discussed so far is the role of the 
atomic SOC parameter a. Even though it is the OAM 
that determines the direction of the angular momentum 
structure, a still plays a crucial role in the Rashba-type 
splitting. Throughout the discussion given so far, we 
considered a large SOC, much larger than the crystal 
field. When a is very small, SAM and OAM do not 
have to be parallel or anti-parallel to each other. This 
means that OAM can point to the same direction for 
two different SAM directions and the state becomes spin 
degenerate (thus the splitting disappears). 

Finally, we argue that, even though ISB at surfaces 
and interfaces (Rashba effect) is considered, our model 
is quite general and has relevance to bulk cases (Dres- 
selhaus effect) [25]. We may expect OAM in the bulk of 
high-Z materials, provided the crystal field can be over- 
come by the strong SOC. In such cases, if there is ISB in 
the bulk like in zinc blend structures, electric dipole inter- 
action can be significant. This may provide microscopic 
mechanism for the large bulk band spin splitting in III-V 
semiconductors [26]. It was also reported that BiTel bulk 
bands have a very large Rashba-type band splitting[27]. 
Since the Bi-layer is sandwiched between Te- and I-layers 
with negatively charged I-layer providing uniform electric 
field to the Bi layer, the situation is similar to the one in 
Fig. 4. Therefore, it will be interesting to check if OAM 



indeed exists in these materials. 
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